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Abstract 

The intestinal microbiota has many beneficial roles for its host. However, the precise mechanisms developed by the 
microbiota to influence the host intestinal cell responses are only partially known. The complexity of the ecosystem and our 
inability to culture most of these micro-organisms have led to the development of molecular approaches such as functional 
metagenomics, i.e. the heterologous expression of a metagenome in order to identify functions. This elegant strategy 
coupled to high throughput screening allowed to identify novel enzymes from different ecosystems where culture methods 
have not yet been adapted to isolate the candidate microorganisms. We have proposed to use this functional metagenomic 
approach in order to model the microbiota's interaction with the host by combining this heterologous expression with 
intestinal reporter cell lines. The addition of the cellular component to this functional metagenomic approach introduced a 
second important source of variability resulting in a novel challenge for high throughput screening. First attempts of high 
throughput screening with various reporter cell-lines showed a high distribution of the response and consequent difficulties 
to reproduce the response, impairing an easy and clear identification of confirmed hits. In this study, we developed a robust 
and reproducible methodology to combine these two biological systems for high throughput application. We optimized 
experimental setups and completed them by appropriate statistical analysis tools allowing the use this innovative approach 
in a high throughput manner and on a broad range of reporter assays. We herewith present a methodology allowing a high 
throughput screening combining two biological systems. Therefore ideal conditions for homogeneity, sensitivity and 
reproducibility of both metagenomic clones as well as reporter cell lines have been identified and validated. We believe that 
this innovative method will allow the identification of new bioactive microbial molecules and, subsequently, will promote 
understanding of host-microbiota interactions. 
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Introduction 

Our body harbors an enormous amount of microorganisms of 
which 90% are found in the distal intestine. This highly dense 
ecosystem has repeatedly been shown to be tighdy associated with 
its host, having a profound and sometimes unexpected impact on 
its health and well-being. If the equilibrated interplay between the 
intestinal microbiota, the host and ingested nutrients is disrupted, 
multiple metabolic, degenerative, inflammatory or infectious 
pathologies can emerge. Over the past years, the development of 
new animal models extended the understanding of the role of 
commensal bacteria to its implication in several physiological 
mechanisms from local effects on the epithelial barrier to systemic 
impacts on immunity and metabolism [1-3]. The intestinal 
microbiota contributes to intestinal homeostasis through direct 
regulation of the development of the intestinal mucosa and 



maturation of the immune system. Therefore the intestinal 
microbiota can be perceived as an integral component of the 
host's physiology [3,4]. 

The complexity of the ecosystem and the difficulty to culture 
most organisms has long prevented in depth functional exploration 
of this neglected organ. The development of molecular methods to 
explore this ecosystem summarized in the term metagenomics 
gave a new boost to gut microbiota research [5-7]. Originally 
based on cloning techniques, the advancement of sequencing 
technologies led to an increasingly profound genomic and genetic 
characterization of this diverse and variable ecosystem. Confirm- 
ing what was known from previous, culture-based studies that the 
intestinal microbiota is composed of a very restricted number of 
phyla (mainly Bacteroidetes and Firmicutes) and that in contrast to 
this uniformity at high phylogenetic levels, there is a considerable 
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inter-individual variability and thus phylogenetic diversity at a 
lower phylogenetic level [8] . 

This new sequence-based techniques to study the intestinal 
microbiota have led to numerous correlation-based studies 
exploring the interrelation of its composition with human health 
and various pathologies like inflammatory bowel diseases (IBD), 
obesity, diabetes and allergies [9,10]. The increased genetic 
accessibility of the intestinal microbiota has further led to a more 
gene-centered understanding of the human microbiota, leading to 
a resurgence in interest for the huge genetic reservoir coded in the 
intestinal microbiome. In order to study the functional potential of 
the intestinal microbiota, a novel approach called functional 
metagenomics has been applied. This approach circumvents the 
challenging cultivation of individual intestinal bacteria by heter- 
ologously expressing the intestinal metagenome in a well-known, 
cultivable host (usually E. coli). Through this elegant approach, 
first studies have identified particular intestinal enzymes such as 
bile salt hydrolases fibrolytic enzymes, xenobiotic metabolizing 
enzymes like P-D-glucuronidases and prebiotic degradation 
pathways confirming the validity of functional metagenomics for 
intestinal bacteria [11-14]. Recent studies by our group have 
proposed to adapt functional metagenomics to study host- 
microbiota interaction in vitro [15,16]. In our first studies, we 
demonstrated that the heterologous expression of the intestinal 
metagenome in E. coli is a valid tool to screen for bioactive 
compounds derived from the intestinal microbiota. Their activity 
patterns have been studied on different eukaryotic cellular models. 
First this combination was validated by studying the effect of 
20725 metagenomic clones on intestinal epithelial cell line 
proliferation [16]. In a second step we established an intestinal 
epithelial reporters assay for NF-kB activation used to screen 2640 
large fragment metagenomic clones for elements activating this 
nuclear receptor of crucial importance for intestinal inflammation 
responses [15]. The observed high variation and low hit rate of 
0.8% when applying stringent criteria indicates that the screening 
of high numbers of metagenomic clones is necessary to identify 
bioactive compounds affecting host gene regulation. This can only 
be achieved in qualitatively exploitable way by a high degree of 
automatisation and adapted analysis techniques as used in high 
throughput screening for drug-discovery [17]. Unlike currently 
reported HTS assays, our screening workflow implicates two main 
independent biological systems: the metagenomic clones, and the 
cell-based assay. Combined they lead to a higher level of 
variability as commonly known from HTS. Cell-based HTS 
usually apply compound libraries to reporter systems [18], while 
HTS of metagenomic clone libraries have been tested for 
enzymatic activity on homogeneous defined substrates [11,12]. 
Screening of metagenomic clones on cellular assays demands the 
control and eventual correction of two screens simultaneously. 
Amplified variability makes screening more chalk-nging and thus 
requires robust and highly reproducible protocols and analysis 
tools. 

In the present study, we describe a robust high throughput 
screening and analyzing method for the screening metagenomic 
clone hbraries on reporter cell lines, adaptable to all the intestinal 
cell-lines tested. This method maximizes metagenomic clone 

expression, minimizes the edge effects and noisiness responsible for 
low confirmation rates observed in previous validation screenings. 
The capacity to screen large numbers (over 10000) of metage- 
nomic clones on different intestinal cell lines enables the 
simultaneous screening of different intestinal metagenomes for a 
variety of mechanisms that genetically influence its host. Mech- 
anisms underlying the complex host-microbiota interaction remain 
poorly understood and are of key importance to design microbiota 



related therapeutic targets. This innovative, robust and adaptable 
method will allow us to identify new bioactive molecular signals 
from the gut microbiota able to influence the host's physiology and 
consequendy the underlying mechanisms at the molecular level. 

Results 

Edge effects with former screening methods 

We first applied the formerly reported screening protocol for 
HTS of metagenomic libraries [15]. Therefore, we developed 
several rc'portcr gene assays (RGA) relevant for the study of host 
cell interactions with the microbiota, including one for ANGPTL4 
and one for PPARy activation. After completing the protocol 
described in the former experimental setup, we noticed a large 
spread of the response distribution in the first eight 96-well 
microtiter plates of the screening. The obtained data showed a 
position-dependent distribution of the signal over the plates on the 
used RGA (Figure lA). The signal \ariati()n was specific to the 
RGA and conserved between plates (data not shown). Figure lA 
shows the measures of the first screening on a PPARy RGA. The 
data points plotted represent the growth of the applied bacteria 
(OD600) versus the reporter gene signal (luciferase activity 
measured in relative light units; RLU) as previously established 
[1,5]. Border and core wells were defined as described in Figure lA 
and will be used for accordingly throughout the manuscript. 
Border and core wells of the screening form visibly distinct groups 
on the plot (Figure lA) indicating a different range of responses for 
the relative parts of the plate. This effect is better visualized when 
border and core data points for all plates of the screening are 
separated and plotted as parallel dot plots side by side (Figure IB). 
Statistical analysis confirmed the significance of this observation 
throughout the tested RGAs (Students t-test: p<0.001 for PPARy 
and NF-kB RGAs, p = 0.0014 for the ANGPTL4 RGA). 
Following these observations, we analyzed and optimized the 
different steps of the screening protocol in order to establish a 
method robust and reproducible enough to be applicable for HTS 
of metagenomic hbraries on different intestinal epithelial RGAs. 

Reduction of cell-related edge effects in RGA 

Reduction of in-homogeneities in the cell monolayers was tested 
using methods previously described on adherent non-intestinal cell 
hnes [19]. Therefore, the cell lines used for the construction of 
RGA assays (HT-29, SW1116 and Caco-2) were homogenously 
seeded and pre-incubated for 0, 0.5 and 1 h at room temperature 
prior to 24 h incubation at 37°C. In accordance to Lundholt el al., 
the in-homogeneities in the cellular monolayer were visibly 
reduced independently of the cell line as observed by light 
microscopy [19]. In the plates without any pre-incubation, the 
wells on the plate extremities showed cell agglomeration at the 
border and non-covered surface areas in their center leading to an 
in-homogenous cellular growth. This effect has been described to 
be caused by convection due to strong temperature gradients at 
the border positions of the plates when moved from room 
temperature to the 37°C incubator [20]. The effect could not be 
reduced by pre-warming the medium (data not shown). The 
homogeneity of cell growth and distribution was quantified using 
crystal violet staining after 24 h incubation at 37°C (Figure 2). 
Agglomeration of the cells at the border of the wells led to an 
increased retention of crystal \ i()l(;t in and between the cells and 
subsequently to a higher signal on the concerned position even 
though the cell did not cover the whole surface of the wells. The 
crystal violet staining technique therefore allowed quantification of 
the edge effect. As illustrated in Figure 2A (left panel), this in- 
homogenous cellular growth effect could be reduced by a pre- 



PLOS ONE I www.plosone.org 



2 



August 2014 I Volume 9 | Issue 8 j el 05598 



An Optimized NTS Method to Study Host-Microbiota Interactions 



A PPARy 
100000 



rolOOOO. 
o 



1000 



laOODDOODODQI 



• Border 
o Core 

• Control 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 
Bacterial Growth (OD 600 nm) 



p=0.0014 




o "bo 




NFkB 
0.8 



0.6 



< CD 



0.2 



0.0 



p<0.0001 




Figure 1. HTS using former protocol. A. Screening of 768 clones (8 
plates) of a metagenomic library for PPARy modulation using HT-29 
reporter gene assay (RGA). Bacterial clone lysates were applied on HT- 
29-PPARy (10% vol/vol) for 24 h and luciferase activity (RLU) was 
monitored. Border, core and control wells are represented in grey, white 
and red dots respectively. B. HT-29-PPARy-luc, SW1 1 1 6-ANGPTL4-luc 
and HT-29-NFkB-SEAP were seeded in 96 well-plates and homoge- 
neously activated with sodium butyrate (2 mM, for PPARy and 
ANGPTL4) or £ coli lysate (for HT-29-NFkB-SEAP) using a pipetting 
robotic workstation. 24 h after activation, luciferase (for PPARy and 
ANGPTL4) and SEAP (for NFkB) activities were monitored. Border and 
core wells are represented in grey and white dots respectively, p-values 
indicated on the figure were obtained by a unpaired t-test of core 
versus border values. 
doi:1 0.1 371/journal.pone.01 05598.g001 

incubation time of 1 h at room temperature, allowing the 
homogeneous settling of the cells in the well before exposing 
them to the convection phenomena when introduced into the 
incubator. The right panel of Figure 2 shows the same data 



represented in box-plots separating core and border wells of the 
plate. The homogenization of the signal with increasing setding 
time at room temperature is clearly visible in both representations 
of Figure 2. Reduction of the edge effect was confirmed by 
statistical testing of the border versus the core wells performing a 
Student's t-test on three independent repetitions (Table 1). The 
transferability of this edge effect reducing technique through 
homogenization of the cell sedimentation was confirmed on 
diflFerent reporter cell-lines (Table 1 and 2). Subsequently the 
responses of reporter systems were assessed in this setup. Figure 2B 
compares two representative PPARy RGA plates homogenously 
activated with sodium butyrate with and without pre-incubation at 
room temperature using a surface plot. The reduction in edge 
effects can clearly be observed in both the surface plot and in the 
box-plot representation of the same data (Figure 2B). Reduction of 
the edge effect was also confirmed by statistical testing of the 
border versus core wells by a Student's t-test on three different 
reporter cell-lines (Table 2). A pre-incubation of 1 h at room 
temperature allowing homogenous settling of the cells, is therefore 
a simple and efficient technique to reduce edge effects that can be 
applied to the used intestinal epithelial RGAs. 

Optimal bacterial growth conditions of metagenomic 
clones for HTS 

Based on the large spread in bacterial growth in the first 
screening described above (Figure lA), we decided to compare the 
formerly reported static growth conditions for the metagenomic 
clones with agitated growth conditions (Figure 3). Four represen- 
tative 96-weU plates with diflFerent metagenomic clones were 
cultured under both, agitated and static culture conditions. Pre- 
cultures were inoculated from the glycerol stock plates and 
cultured for 24 h at 37°C. The subsequent cultures were 
monitored by regxilar measurements of the optical density and 
are represented for one representative 96 well plate in figure 3B. 
The growth curves for each individual clone are represented by 
connected grey dots. The plates means for each time-point are 
represented as red dots. The growth curves of agitated culture 
show the canonical exponential growth phase followed by a 
stationary phase. We found that growth was more homogenous in 
agitated as compared to the static culture conditions, resulting in a 
narrower range of distribution of optical density (OD600) and 
consequently bacterial cell density in the plate (figure 3B). In 
addition, static growth conditions show no exponential growth 
phase, which might impact expression of tire metagenomic clones 
insert (Figure 3B). Figure 3A shows box-plots of border and core 
data points for both agitated and static culture conditions. This 
representation further reveals that static growth conditions exhibit 
strong edge effects (p-value calculated by Student's t-test is 
p = 7.8x 10"''*) while agitated culture shows no difference in 
growth distribution between core and border positions (p- 
value = 0.1 1). Agitated growth condition further highlight outliers 
due to impaired growth that are covered by the distribution under 
static growth conditions. These data clearly show that in order to 
achieve a better, homogenous bacterial growth associated with the 
elimination of edge effects, the use of agitated culture conditions is 
more appropriate for HTS since it guaranties the comparability of 
the single tested metagenomic clones. 

Probability for heterologous expression can be increased and 
insert expression levels can be homogenized through optimization 
and homogenization of bacterial growth in all wells. Further most 
cell lines responses are biased, though in different extend by the 
growth level or the bacterial host of the metagenomic library. 
Therefore homogenization of metagenomic clones growth also 
means a homogenization of background noise. 
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Table 1. Cellular growth edge effect reduction. 







Mean p-values + SEM 


T = 0 


T = 1 h 


SW1116 


0.014+/-0.012 


0.202±0.089 


HT-29 


0.023 ±0.023 


0.063 ±0.062 



Cell-lines were processes as in Figure 2 using crystal violet. Mean p-values -t-/— 
SEM comparing borders versus core for different parental cell-lines at 0 or 1 h 
room temperature incubation are reported. 
doi:l 0.1 371/journal.pone.0105598.t001 



Homogenization of metagenomic insert expression and 
activity for HTS 

To screen metagenomic clone libraries for bioactive compounds 
we need to guaranty the expression of the insert and the 
avaUabihty of the bioactive compound. First can be done by 
defining the optimal expression conditions while for the second a 
lysis protocol, that does not impact bioactivity has to be 
estabhshed. In order to estabhsh these protocols that guaranty 
homogenous growth, metagenomic insert expression and avail- 
ability, we used a metagenomic clone named 43C8 exhibiting 
partially secreted alkaline phosphatase activity. Growth was 
monitored by the measurement of the OD600 while metagenomic 
insert expression was quantified using a colorimetric assay to 
monitor alkaline phosphatase activity (QUANTI-Blue reagent, 
OD655). The 43C8 metagenomic clone was grown in 96-well 
plates under agitated and static conditions. OD600 was monitored 
and samples were frozen at every time point prior to subsequent 
measure of the phosphatase activity. Data shown in figure 4A 
confirm suboptimal growth under static conditions as opposed to 
exponential growth under agitated growth conditions (700 RPM). 
For the metagenomic clone 43C8, it further confirms the 
hypothesis that optimal growth supports and homogenizes the 
expression of metagenomic inserts. Agitated culture shows 
exponential augmentation of phosphatase activity over time while 
static growth leads to a slower augmentation of activity and a 
higher spread within the four independent wells used for the 
experiment. Figure 4A further suggests that metagenomic insert 
expression can be delayed under unfavorable growth conditions. 
We therefore suggest a culture time of at least 1 6 h covering all 
bacterial growth phases to reach full expression of the metage- 
nomic insert circumventing potential growth phase dependent 
expression. 

Metagenomic clones used for screening carry a 40 Kb insert of 
heterologous DNA. Even though whole operons have been 
described in metagenomic clones, the secretion of eventual 
activating compounds cannot be presumed. We therefore decided 
to lyse the metagenomic clones in order to solubilize intracellular 
compounds. 



In order to conserve integrity of eventual bioactive compounds 
lyses using freeze thaw cycles was chosen. As represented in the 
parallel dot plot of entire 96-well microtiter plates in Figure 4B, 
lysis not only homogenizes the activity of the single clones within 
the plate but also renders the intracellular compounds accessible 
thereby increasing the signal of the 43C8 clone. Even though the 
repetition of freeze/ thaw cycles further homogenizes metagenomic 
clone activity, two freeze thaw cycles seem to be the most adequate 
lysis for HTS application showing a robust result with manageable 
effort. 

Bacterial and cellular growth optimal protocols reduce 
the background noise and increase the reporter assay 
sensitivity 

The newly estabhshed protocols were combined and tested for 
their effect on reporter assay sensitivity. Intestinal epithelial cell 
lines display differential amounts of activation or modification 
upon application of lysates of £. coli such as used for the functional 
metagenomic screening. RGA used to screen host-microbiota 
interactions, such as NF-kB, Angptl4 and PPARy are known to be 
activated by bacterial compounds other than coming from 
metagenomic insert. Therefore the background noise for the 
screening is a combination of RGA noise and bacterial growth 
homogeneity. Thus we tested the sensitivity of screen when 
applying the previously established modifications. We randomly 
applied three serial dilutions of the PPARy activating component 
sodium butyrate (NaBut; 500 |xM, 1 60 |aM and the sub-activating 
concentration 80 |J,M) equally to border and core positions of a 
RGA plate using the previously mentioned optimization tech- 
niques. The HT-29 PPARy reporter cell line was seeded and pre- 
incubated for 1 h at room temperature before 24 h incubation at 
37°C. Prior to NaBut stimulation, all wells were homogenously 
activated with a filtered lysate of the E. coli strain culture used for 
the metagenomic library, namely EPI300-Cont, to simulate 
screening conditions. Figure 5 represents the PPARy response to 
the three NaBut concentrations from three independent plates. 
The position of the individually activated wells is indicated next to 
the plotted points. Wells with activator concentrations of 500 |.lM 
(EC50, Figwe SI) and 160 [iM NaBut were clearly identifiable. 
Interestingly, even with a high background signal related to the 
control E. coli EPI300-Cont stimulation, we were able to identify 
activation with a concentration of NaBut at its detection limit 
(80 |J,M) for most of the wells. Assuming a normal distribution of 
data-points the threshold for hit detection is set at 3 standard 
deviations from the mean (mean-l-3 SD). A normalization versus 
the mean as illustrated in Figure 5 allows the detection of clear 
signals. The cutoff of 3 standard deviations however generates falls 
positives and negatives even at this small scale and with known 
activators. We therefore tested methods to minimize position 



Table 2. Cellular growth-dependent edge effect reduction applied to three different reporter cell-lines. 





Mean p-values +/— SEM 


T = 0 


T = 1 h 


Activator 


HT-29 NF-kB-SEAP 


0.002±0.002 


0.248±0.147 


E coli lysate 


SW1116 ANGPTL4-IUC 


0.000079 ±0.00008 


0.202 ±0.068 


Sodium Butyrate {2 mM) 


HT-29 PPARy-luc 


0.023±0.015 


0.239±0.108 


Sodium Butyrate {2 mM) 



Cell-lines were processes as in Figure 3. The table represents the mean of three p-values +/— SEM for comparing borders versus core for the different tested reporter 
cell-lines. 

doi:l 0.1 371/journal.pone.Ol 05598.t002 
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Figure 2. Edge effect reduction on cellular growth applied to 
the HT-29-PPARy reporter cell-line. A Parental HT-29 cells were 
homogeneously seeded using a pipetting robotic workstation and pre- 
incubated at room temperature for the indicated time (0, 0.5 or 1 h) 
prior to 37'C incubation for 24 h. Cellular monolayer homogeneity was 



monitored using cell staining with crystal violet and quantified by 
absorption measurement at 595 nm. Left panel is a surface plot of a 
representative set of plates at increasing times of room temperature 
pre-incubation. Boxplot representation (right panels) of the OD 595 
values for the border compared to the core of the respective 
representative plates. Mean p-values for different cell-lines tested are 
summarized in Table 1. B HT-29-PPARy reporter cells were homoge- 
neously seeded using a pipetting robotic workstation and pre- 
incubated at room temperature for the indicated time (0 or 1 hour) 
then at 37 C for 24 h prior to sodium butyrate (2 mM) activation. 
Luciferase activity (RLU) was quantified after 24 h of activation. Graph 
(left panels) shows two representative plates with different pre- 
incubation times at room temperature. Boxplots (right panels) 
represent the border and core values of the surface plots. Mean p- 
values for different cell-lines tested are summarized in Table 2. 
doi:10.1 371/journal.pone.01 05598.g002 

effects and clearly eliminate false positives and negatives in this 
controlled setup. 

Normalization as an additional tool to increase the 
robustness of functional metagenomic HTS 

In order to increase robustness and sensitivity of our assay, 
different data treatment strategies were tested. Since protocol 
optimization has reduced edge effects dramatically, but does not 
eliminate position effects completely in all the reporter ceU-lines 
used (Figures 1-5 and data not shown), we tested different data 
analysis protocols for their capacity to reduce positional effects 
found in all biological screenings and so increase the sensitivity and 
robustness of our assays. After homogenization of metagenomic 
clone culture conditions and thus elimination of expression bias 
due to growth, our readout can be considered to be one 
dimensional, as opposed to the two dimensional view defining 
each data point with the RGA signal and the metagenomic clones 
culture OD600. We consequently decided to present the data as 
ranked spots in order to better visualize cut-ofiFs at both extremes 
(Figure 6). Raw data from the reported first screening on HT-29- 
PPARy were taken for analysis (Figure lA). Figure 6 A shows the 
ranked untreated data where red dots represent data points from 
the core wells and green points represent data points from the 
border wells of the 96-well plates. Superposition of green and red 
points results in dark green points, indicating the part of the curve 
where border and core data point overlap. The separation of core 
and border data points showed previously in Figure lA is also 
clearly visible in this representation (Figure 6A) where both 
extremes are dominated by one color. The frequentiy used 
normalization method Z-score facilitates inter plate comparison, 
but does not eliminate position effects of the plates as the clear 
separation of border and core in Figure 6B illustrates. To consider 
the observed position effects we used a more appropriate 
normalization method less susceptible to outliers and position 
effects namely the B-score. Roughly in order to take into account 
column and row effects the B-score represents the ratio of an 
position adjusted raw value to a measure of variability in the plate 
named Median absolute deviation (MAD) [18,21,22]. The 
resulting normalized data are shown in Figure 6C. B-score 
normalization reduces the positional effects and simultaneously 
accentuates the tails of the distribution indicating that B-score 
normahzation is an appropriate tool to address the position effects 
in our system. We then applied the B-score normalization to the 
dataset represented in Figure 5 where randomly-picked wells of 
the core and border of the plates were stimulated with RGA 
activating NaBut down to its previously defined activation limit 
(Figure SI). This dataset generated using the optimized protocol 
had a lower spread of its data point distribution (Figure 7). The 
reduction of variation visibly improved the detection of activators. 
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Figure 3. Metagenomic clones growth optimization. Four random plates of EPI300 metagenomic clones were inoculated from their -80°C 
glycerol stock into 96 well-plates (10% v/v) using a pipetting robotic workstation. After a 24 h pre-culture, they were re-inoculated and cultured 
overnight (1 0% v/v) at 37 C. These steps were performed either in static or agitated (700 rpm) culture conditions. Bacterial growth was monitored by 
optical density (OD 600 nm). A. Border and core wells of bacterial cultures in agitated (left panel) or static culture-conditions (right panel) for thus 
representative set of 4 plates, p-values indicated on the figure are the result of a unpaired t-test of the core versus border bacterial culture ODs. B. 
Growth curves of one representative plate of metagenomic clones (OD 600 nm) for agitated and static culture conditions. Red dots: mean OD600, 
joined grey dots: growth curve of single metagenomic clones. 
doi:1 0.1 371 /journal.pone.01 05598.g003 



Ranked raw data shown in Figure 7A already allows the detection 
of strong activators quite robustly at a small scale, tough low 
activators are partially covered by the variation of the background 
noise. Z-score normalization allows an inter-plate normalization 
but does not increase sensitivity of the assay since the Z-score is 
sensitive to outliers; in this case, the activated wells (figure 7B). 
The activated wells increase variability and thus the recommended 
cut-off of mean ±3 SD used to define hits covers a large part of 
activators. In contrast, B-score normalization is insensitive to 
position effects. Applied to the same dataset, Figure 7C shows that 
B-score normalization cannot only smoothen the distribution 
through its elimination of position effects. It also accentuates the 
outliers, and allows the low concentration of activators to be 
clearly identified dirough a clear-cut separation from the 
background noise. For our screening setup, B-score normalization, 
therefore, avoids the loss of hits as false negatives due to 
background noise and simplifies the identification of hits through 
its accentuation of outliers. The combination of an experimental 
reduction in variation completed with statistical elimination of 
position effects resulted therefore in a highly sensitive assay. 



meeting all requirements for a High Throughput Screening assays 
on host-microbiota interactions. 

Screening using the normalization protocol for 
identification of metagenomic clones responsible for 
bioactive effects 

The established protocol was used to screen 92 microtiter-plates 
with a total of 8832 clones (Figure 8). Analysis of Z-score 
normalized and B-score normalized data confirms our previous 
observations that, even with an optimized protocol reducing 
position effects in all steps of the screen, there remains a position 
bias. Figure 8A where data points on the border are plotted in 
green and wells of the core are plotted in red illustrates the position 
bias Z-score normalization propagates. This makes it impossible to 
identify clones using the common criteria of intra-plate Z-score 
normalization. 

B-score in contrast showed a much narrower distribution where 
inactivated wells are concentrated to the mean (Figure 8B). This 
distribution with its accentuation of oufliers allows a clear 
separation of background-noise and positives on the ranked plot 
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Figure 4. Bacterial growth, insert-dependent activation and lyses process optimization. The alkaline phosphatase-expressing clone 43C8 
were cultured under agitated conditions, harvested at different times of growth, and subsequently lysed by congelation cycles (>1 h at 80'C). After 
filtration on 0.45 |im filters phosphatase alkaline substrate (QUANTI-Blue) was activated with bacterial lysates from different culture time-points. After 
24 h, alkaline phosphatase activity was monitored (SEAP activity at OD655). A. represents the growth curve of the 43C8 (OD600) and the phosphatase 
alkaline activity of the respective lysates (AP activity at OD655). B. represents the alkaline phosphatase activity (OD655) under different lyses cycles 
(Control represents £ coU bearing an inactive fosmid with a human 40 kb-insert). 
doi:1 0.1 371/journal.pone.01 05598.g004 



and thus a better identification of iiits. Therefore the established 
method combined with B-score normalization confirms to be 
suitable for HTS of metagenomic libraries on RGAs. 

The strategy for functional metagenomic assays aims to the 
identification of the bacterial genes involved in the metagenomic 
clone's activity on the selected reporter assay. Once hits are 
identified in the screening downstream processing of confirmation 
screening and identification of bioactive compounds wiU be 
performed in a small scale. The presented method is thus a 
powerful tool to identify bioactive compounds of the intestinal 
microbiota. Downstream processing however has to be adapted to 
the nature of the identified activator and is thus not subject of this 
study [15]. 



Discussion 

All increasing number of studies emphasize the diverse roles of 
the intestinal microbiota. These roles are often associated with 
beneficial physiological effects for the host, reflecting a symbiotic 
and beneficial crosstalk between cells of the intestinal epithelium 
and the resident microbiota. However, the cellular and molecular 
mechanisms developed by the microbiota in order to influence the 
host's intestinal cell responses remain poorly understood. 

An emerging and innovative field called functional metage- 
nomics, has led to the discovery of a variety of novel genes and 
their associated functions. This approach was first applied to the 
intestinal microbiota to identify enzymes such as bile salt 
hydrolases, fibrolytic and xenobiotic metabolizing enzymes [1 1— 
1 3] . This method has proven to be effective since the well-known, 
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Figure 5. Detection limit on optimized reporter gene assay. The HT-29-PPARy cell-line was homogeneously seeded using a pipetting robotic 
workstation and incubated at room temperature for 1 h prior to 37 C incubation. After 24 h of culture, all wells were activated homogenously with 
EPI300-Cont lysates for 24 h at 37'C. In randomly chosen wells (50% border and 50% core), three different concentrations of PPARy activating NaBut 
(at 500 liM, 160 |ilVl and the sub-activating concentration 80 |iM) were randomly added in the 96 well-plates. The graph represents normalized values 
from 3 independent plates. 
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cultivable E. coli, used as host, is able to express up to 40% of the 
functional potential from randomly cloned environmental DNA. 
Among others, new enzymes from anaerobic and Gram-positive 
bacteria have been identified by heterologous expression of 
candidate metagenomes in the easy cultivable E. coli [23]. We 
previously proposed the use of a functional metagenomic 
approach to model the host-microbiota interaction in vitro 
[13,15]. Importantly, previous functional studies showed that hit 
rates in metagenomic screenings are expected to be in the range of 
less than 0.01%, therefore requiring high throughput screening in 
order to identify molecules of interest. Previous metagenomic high 
throughput screenings (HTS) were performed to identify enzy- 
matic activities on well-defmed substrates. Contrary to these first 
metagenomic screenings, the study of host cell-microbiota 
interactions described here has to take into account two 
independent biological systems: the metagenomic clones, and the 
cell-based assay, representing two independent sources of 
variability and noise for the screen. This was observed by using 
the former experimental setup on different reporter cell-lines 
(Figure 1 and data not shown) where we noticed a position- 
dependent and highly spread distribution of the response. These 
sources of variations and consequent difficulties to reproduce the 
observations are not compatible with HTS since they impair the 
ability to robustly identify hits, leading to high numbers of false 
positive and negatives. We therefore developed a robust and 
highly reproducible protocol and defined appropriate analysis 
tools in order to be able to exploit this promising new approach on 
a broad range of reporter-cell line assays. 

In the present study, we established the conditions necessary to 
perform a robust, reproducible high throughput screening assay 
combining metagenomic clones and cellular models. This allows 
us to identify bioactive motives from the human intestinal 
microbiota on different intestinal cell-line models, and hence 
screen the intestinal microbiome for the single and eventually rare 



mechanisms implicated in host-microbiota cross-talk. In this 
optimal setup, we minimized the edge effects, variation and 
noisiness of the assay, we had previously identified as the source for 
high insecurity level for hit definition and subsequendy low 
confirmation rates. The optimal methods implies the use of a 
robotic pipetting workstation, a 1 6 h agitated bacterial culture and 
a 1 h room temperature settling time for the reporter cell-line. The 
increased control of the screening process is completed by 
appropriate data treatment, namely the B-score normalization, 
which showed to be an appropriate method to normalize increased 
the robustness and sensitivity to a level compatible with HTS. 
These adaptations provided clear identification of minimal 
activator concentrations with a threshold of as little as one 
standard deviation, which demonstrates that HTS can be 
performed on combinations of two biological systems like the 
combination of metagenomic clones with cell models to identify 
interactions. We believe that the presented approach as summa- 
rized in Figure 9 can be used to screen different gene pools 
reaching from complex ecosystems microbiomes to genomic 
approaches for single microbes. 

The understanding of the complex host-microbiota relationship 
and the possibility to modulate these interactions is of key 
importance for human health. Our knowledge on the various 
contributions of the microbiota to health is still in its infancy and 
the underlying cellular and molecular mechanisms of its interplay 
with the host intestinal cells remain poorly understood. We believe 
that this innovative and adaptable method will allow the 
identification of new bioactive molecules from the intestinal 
microbiota able to influence the host's physiology on various levels 
such as immunity, metabolism or tissue differentiation. Using this 
optimized method, we have already identified over 30 clones of 
interest which are being fully characterized using the additional 
experimental steps previously described [15]. Expanding this new 
form of HTS to new targets and ecosystems is a challenge that wiU 
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Figure 6. Different normalizations of assay plates. Different 
mathematical approaches for normalization were applied to the dataset 
from the initial HT-29-PPARy screening. A shows all data-points as dots 



ordered by activity-rank whereby the red dots are points on the plate's 
cores and the green points are well on the borders. Dark green areas 
show overlapping of red and green dots. B shows ranked Z-score 
normalization of the results, (red dots are points on the cores and the 
green points are well on the borders. Dark green areas show 
overlapping of red and green dots. C. shows the ranked B-score values 
for the same data set. (red dots are points on the cores and the green 
points are well on the borders. Dark green areas show overlapping of 
red and green dots. 
doi:10.1 371/journal.pone.01 05598.g006 

build on this first experience. Identifications of the bacterial actors 
of such interplay will lead to the molecular characterization of the 
mechanisms underlying the microbiota's properties. This ap- 
proach wUl allow the identification of molecular signals from yet 
totally undescribed micro-organisms as well as novel mechanisms 
of the host responses. Considering the importance of host- 
microbiota interactions for human health, the access to yet 
unexplored resources will lead to the discovery of new actors 
regulating this cross-talk and therefore new potential therapeutic 
targets. 

Materials and Methods 

Cell Culture and Reagents 

HT-29 (ATCC HTB-38) and SW1116 (ATCC CCL233) cells 
were from ATCC and maintained in RPMI 1640 medium in a 5% 
C02 atmosphere at 37°C. Both media were complemented with 
10% fetal calf serum, 10 mM sodium pyruvate, 10 niM nones- 
sential amino acids, 2 mM glutamine, penicUlin (50 U/ml), and 
streptomycin (.50 U/ml) (all from Lonza). 

Plasmid Constructs 

The ANGPTL4 promoter sequence was amplified from human 
genomic DNA using the HiFi DNA polymerase (Fermentas) 
according to the manufacturer's instructions. Briefly, primers used 
for ANGPTL4 PCR-amplification were fwiTTTTTTGCTA- 
GCCTCAGGACATTAAAGACCCTGGCGGTAGAG and rw: 
TTTTTTGGTACCCCTCTTAGGTAGCCTGGGAGCGGG- 
GATTCG. The ANGPTL4 promoter construct (consisting of the 
1 .8 Kb upstream transcriptional start) was cloned into a modified 
pcDNA3.1/Zeo (harboring luciferase and without any CMV 
promoter sequence, kind gift from A. Cultrone, INRA, France) 
and verified for authenticity by sequencing [24]. PPARy reporter 
construct pJ3-TK-Luc was from M. Chamaillard, (INSERM Lille, 
France) pTK-Hygro was from Invivogen. 

Transfection and stable cell-lines selection 

Transfections were performed using TFX50 (Promega), accord- 
ing to manufacturer's recommendations. Stable reporter cell lines 
for PPARy were cotransfected with pTK-Hygro (Clontech) and 
the ANGPTL4 promoter sequence, were selected by selective 
pressure of the appropriate antibiotic (600 |J,g/ml Hygromicin, 
50 Hg/ml Zeocin). The NF-kB reporter clone HT-29/kb-seap-25 
has been described previously [15]. 

Reporter cell-lines activation and analysis 

All experiments were performed using a robotic pipetting 
workstation (Hamilton Robotics Star Line) for high seeding, 
stimulation and revelation homogeneity. For each experiment, 
cells were seeded at 25000 cells/ well in 96 well-plates and incubated 
24 h prior to activation. Homogeneous cellular activations were 
performed using sodium butyrate (NaBut, 2 mM, Sigma) for 
PPARy and ANGPTL4 reporter systems or 0,45 jim-filtered, 
sonicated E. coli (Epi300 from Epicentre, 10% vol/vol) for 
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Figure 7. Detection limit of the screening protocol. Statistical 
normalization applied to optimized cellular and bacterial culture 
conditions. As in Figure 5, HT-29-PPARy cell-line was homogeneously 
seeded using a pipetting robotic workstation and incubated at room 
temperature for 1 h prior to incubation. After 24 h homogenous EPI300- 
Cont lysates were applied to the cell assay (1/10 v/v) for 24 h at 37 c. In 
randomly chosen wells (50% border and 50% core), three different 
concentrations of PPARy activating NaBut (at 500 |jM, 160 |ilVl and the 
sub-activating concentration 80 |ilVl) were randomly added in the 96 
well-plates. The resulting measured luciferase activity (expressed in 
relative light units, RLU) is represented as activity-ranked dots. Luciferase 
activation for EPI300-Cont, EPI300-Cont4-500 |ilVI NaBut, EPI300-Cont4- 
160 liM NaBut and EPI300-Cont4-80 |ilVI NaBut are represented in black, 
red, green and blue respectively. Grey bars represent mean activation 
(Mean) and mean activation 4-/— standard deviation (SD). 
doi:10.1 371/journal.pone.01 05598.g007 

NFkB-SEAP system. Cells were then incubated 24 h prior 
quantification of the reporter gene expression (alkaline phosphatase 
(SEAP) or luciferase). SEAP (for NFkB system) secretion or 
intracellular luciferase expression (for both ANGPTL4 and PPARy 
systems) were revealed using the QUANTI-Blue (Invivogen) or 
One-Glo (Promega) reagent respectively, according to the manu- 
facturers' instructions. AU measurements were performed using a 
microplate reader (Infinite 200, Tecan) and quantified at 655 nm 
OD for SEAP and by luminescence for luciferase. 

Crystal violet staining 

Cellular growth was monitored using crystal violet absorption 
property at OD 595 nm as previously published [25] . Briefly, HT- 
29 or SW1116 cell-lines were homogeneously seeded using a 
pipetting robotic workstation and incubated at room temperature 
for the indicated time (0-1 h) prior to 37°C incubation for 24 h. 
Then, cells were fixed in methanol at — 20°C for 15 min prior 
incubation with 0. 1 % of crystal violet in deionized water for 30 min. 
Subsequentiy, cells were extensively washed in water. Crystal violet 
retained in the cell layer was solubUized in 100 |.ll 0.2% Triton X- 
100 (Sigma) and quantified through absorption measurement at 
595 nm using a microplate reader (Infinite 200, Tecan). 

Metagenomic clones 

Randomly selected plates from the metagenomic library issued 
from the intestinal microbiota of obese and lean patients 
(constructed using CopyControl Fosmid Library Production Kit 
by Genoscope, France) were used in this study. The library consists 
of E. coli EPI300T1R clones carrying each one fosmid with a 
metagenomic inserts of approximately 40 kb. E. coli bearing an 
inactive fosmid with a human 40 kb-insert (ligation control 
reaction for CopyControl Fosmid Library Production Kit), 
referred to as EPI300-Cont, was used as control. Bacteria were 
grown under static or agitated conditions in 96 weU-plates in LB 
broth (Lennox L broth base, Invitrogen) at 37°C. Bacterial growth 
was monitored by absorbance measurement at 600 nm, using a 
microplate reader (Infinite 200, Tecan). We used a metagenomic 
cloneshowing partially-secreted alkaline phosphatase activity, 
namely 43C8, as a read-out for bacterial lysate activity. Bacterial 
lysates were obtained through freezing at — 80°C and subsequent 
thawing at room temperature. Bacterial lysates were filtered by 
centrifugation through a 0.45 |am microplate filter (Corning) prior 
to use on reporter cell-lines. 

Statistical analysis 

Graphic representations and statistical analysis were performed 
using the Graphpad Prism (GraphPad Software, San Diego 
California USA) or R statistics softwares (www.r-project.org/). 
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Figure 8. PPARy HTS screening. Comparison of different data analysis methods of screening data from the screening of 92x96 metagenomic 
clones on a PPARy RGA. Figure A represents Z-score normalized data B represents B-score normalized data. The top left picture plots the each 
readout in function of the growth of the metagenomic clone. The bottom left picture represents the distribution of the readouts over the whole 
dataset and the large picture on the right part represents all data points ranked by readout. The coloration of red dots for the borders and empty 
green dots for the core allow to visualize eventual position effects. 
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Functional Metagenomic Library 



•Optimal growth parameters 
•Intra- and interplate growth homogeneity 
•Metagenomic insert expression 
maximization 



•Agitated growth for 16h in microtiter 
plates 

•Lysis through 2 freeze thaw cycles 
•Filtering through 0.45 pm filter 
•Application of 5-10 pi per well 




Intestinal Reporter Gene Assay (RGA) 



•Adaptation to different reporter genes 
(luciferase and SEAP) 
•Protocol for reduced position effect on 
different intestinal cell lines 
•Validation on HTS scale 

•Establish and characterize RGA for target 
gene in intestinal epithelial cell line 

•Seed in appropriate dilution, settle for 1h 
at room temperature 

•Incubate for 24h at 37°C prior to activatior 



Screening and hit identification 

Validation of analysis of new data type 
Automatization 

Stimulate RGA for the required time 
Reveal readout and retrieve data 
Normalize data using B-score 
normalization 

Identify hits on ranked B-score plot 



Adaptation of the Protocol 
Final HTS Protocol 



Figure 9. Schematic overviewed of challenges addressed and the resulting suggested protocol. 

doi:10.1371/journal.pone.0105598.g009 



Comparisons of distributions were performed using an unpaired t- 
test with 95% confidence intervals, p Values smaller than 0.05 (p^ 
0.05) were considered significant. Z-scores and B-scores were 
calculated as described by Malo et al. using R statistics software 
[18]. Median polish was performed plate-wise based on the 
Tukeys method [26] using the R command medpohsh(x). 
Residuals were used for further calculation of the B-score 
(medpolish(x)$residuals). 

Supporting Information 

Figure SI Dose response curve of sodium butyrate on 
the PPARy RGA used in this study. 

(EPS) 



Acl(nowledgments 

Wc thank the MetaFun platform within MetagGenopolis (INRA, Jouy-en- 
Josas) for the use of the robotic pipetting workstation. We also thank M. 
Chamaillard, A. Cultrone and O. Lakdhari for useful constructs and 
adxdces. We thank Laure Laval for proxdding the control E. coli strain 
(EPLSOO-Contr). 

Author Contributions 

Conceived and designed the experiments: TDWJD HMB NL. Performed 
the experiments: TDW MN NL FL. Analyzed the data: TDW HMB NL. 
Contributed to the writing of the manuscript: TDW^ HMB NL. 



References 

1. Sekirov I, RusscU SL, Antunes LCM, Finlay BB (2010) Gut Microbiota in 
Health and Disease. Physiol Rev: 8,59-904. doi:10.11,52/physrcv.00045.2009. 

2. Hooper L V, Midtvedt T, Gordon JI (2002) How host-microbial interactions 
shape the nutrient environment of the mammalian iutestiue. Annu Rev Nutr 22: 
283-307. AvaQablc: http://wvvvv.ncbi.nlm.nih.gov/pubmed/ 120.').5347. Ac- 
cessed 2011 Jun 10. 

3. Geuking MB, KoUcr Y, Rupp S, McCoy KD (2014) The interplay between the 
gut microbiota and the immune system. Gut Microbes 5. Available: http:// 
www.ncbi.nlm.nih.gov/pubmed/24922519. Accessed 2014 Jul 16. 

4. Ivanov II, Honda K (201 2) Intestinal commensal microbes as immune modulators. 
Cell Host Microbe 12: 496-508. Available: http://www.pubmcdcentral.nih.gov/ 
articlerender.fcgi?aitid = 351 6493&tool = pmccntrcz&rendertvqx' = abstract. Accessed 
2014JuI 15. 



Suau A, Bonnet R, Sutren M, GodonJJ, Gibson GR, et al. (1999) Direct analysis 
of genes encoding 16S rRNA from complex communities reveals many novel 
molecular species within the human gut. Appl Environ Microbiol 65: 4799— 
4807. Available: http:/ /www. pubmedcentral.nih.gov/articlerender. 
fcgi?artid = 91647&tool = pmcentrez&rendertype = abstract. 
Hayashi H, Sakamoto M, Benno Y (2002) Phylogenetic analysis of the human 
gut microbiota using 16S rDNA clone libraries and strictly anaerobic culture- 
based methods. Microbiol Immunol 46: 535—548. Available: http://www.ncbi. 
nlm.nili.gov/pubmcd/12363017. Accessed 2011 Apr 21. 

Handelsman J (2004) Metagenomics: Application of Genomics to Uncultured 
Microorganisms. Microbiol Mol Biol Rev 68: 669-685. doi:10.1 128/MBR. 
68.4.669. 



PLCS ONE I www.plosone.org 



12 



August 2014 | Volume 9 | Issue 8 | e105598 



An Optimized HTS Method to Study Host-Microbiota Interactions 



8. QinJ, Li R, RacsJ, Arumugam M, Burgdorf' KS, ct al. (2010) A human gut 
microbial gene catalogue established by mctagenomic sequencing. Nature 464: 
59-65. Available: http://www.ncbi.nlm.nih.gov/pubmed/2O2036O3. 

9. Lepage P, Leclerc MC, Joossens M, Mondot S, Blottiere HM, et al. (2012) A 
metagenomic insight into our gut's microbiome. Gut. Available: htip:/ / gut.bmj. 
com/cgi/doi/10.1136/gutjnl-2011-301805. Accessed 2012 Apr 24. 

10. Blottiere HM, dc Vos WM, Ehrlich SD, Dorc J (2013) Human intestinal 
metagcnomics: state of the art and future. Ciwr Opin Microbiol 16: 232-239. 
Available: http://w\vw.nrbi.nlm.nih.go\'/ pubmed/23870802. Accessed 2013 
Aug 8. 

1 1 . Jones B V, Begley M, Hill C, Gahan CGM, Marches! JR (2008) Functional and 
comparative metagenomic analysis of bile salt hydrolase activity in the human 
gut microbiome. Proc Natl Acad Sci U S A 105: 13580-13585. Available: 
http:/ /www. pnas.org/ cgi/ content/ abstract/ 105/36/1 3580. Accessed 201 1 
Aug 20. 

12. Tasse L, Bcrcovici J, Pizzut-Serin S, Robe P, Tap J, ct al. (2010) Functional 
metagcnomics to mine the human gut microbiome for dietary liber catabolic 
enzymes, (ienome Res 20: 1605—1612. Available: http://\vww.ncbi.nlm.nih. 
gov/pubmcd/20841432. Accessed 2010 Sep 18. 

13. Gloux K, Berteau O, El Oumami H, Beguet F, Leclerc M, et al. (2011) A 
metagenomic P-glucuronidase uncovers a core adaptive fimction of the human 
intestinal microbiome. Proc Natl Acad Sci U S A 108 Suppl : 4539-4546. Available: 
http:/ / www.pubmedcentral.nih.gov/ articlercnder.fcgi?artid = 3063586&tool = pmcentrez 
&xendertype = abstract. Accessed 201 1 Jul 29. 

14. Cccchini DA, Laville E, Laguerre S, Robe P, Leclerc M, ct al. (2013) Functional 
metagcnomics reveals novel pathways of prebiotic breakdown by human gut 
bacteria. PLoS One 8: e72766. Available: http://www.pubmedcentral.nih.gov/ 
articlerender.fcgi?artid = 3774763&tool = pmcentrez&xendertype = abstract. Ac- 
cessed 2014 Apr 7. 

15. Lakhdari O, Cultrone A, Tap J, Gloux K, Bernard F, et aL (2010) Functional 
Metagenomics: A High Throughput Screening Method to Dec^her Microbiota- 
Diiven NF-kB Modulation in the Human Gut. PLoS One 5: el3092. Available: 
http://www.pubmedcentral.nih.gov/ article re ndcr.fcgi?artid = 2948039&tool = 
pmcentrez&rendertype — abstract. Accessed 2010 Oct 1 . 

16. Gloux K, Leclerc M, Iliozer H, L'Haridon R, Mamchanh C, et al. (2007) 
Development of high-throughput phenotyping of metagenomic clones from the 



human gut microbiome for modulation of cukaryotie cell growth. Appl Environ 
Microbiol 73: 3734-3737. Available: http://www.pubmedcentral.nih.gov/ 
articlerender.fcgi?artid = 1932692&tool = pmcentrez&rendertype = abstract. Ac- 
cessed 2011 Jun 22. 

17. Macarron R, Banks MN, Bojanic D, Bums DJ, Cirovic D, et al. (201 1) Impact of 
high-throughput screening in biomedical research. Nat Rev Drug Discov 10: 
188-195. Available: http://www.ncbi.nlm.nih.gov/pubmed/21358738. 

18. Malo N, HanleyJA, Cerquozzi S, Pelletier J, Nadon R (2006) Statistical practice 
in high-throughput screening data analysis. Nat Biotechnol 24: 167-175. 
Available: http://dx.doi.org/io.1038/nbtll86. Accessed 2011 Jul 18. 

19. Lundholt BK, Seudder KM, Pagliaro L (2003) A simple technique for reducing 
edge effect in cell-based assays. J Biomol StMXTu 8: 566 570. Available: http:// 
w\\w.nebi.nlm.nih.gov/pubmcd/14r)6778 1. Accessed 2010 Jul 23. 

20. Esscr P (2010) 'I'hermo Scientific Nunc Edge Plate 'I'echnical Information. 

21. Makarenkov V, Zentilli P, Kcvorkov D, Gagarin A, Malo N, et al. (2007) An 
efficient method for the detection and elimination of systematic error in high- 
throughput screening. Bioinformatics 23: 1648-1657. Available: http://www. 
ncbi.nlin.nih.gov/pubmed/17463024. Accessed 2011 Mar 28. 

22. Brideau C, Gunter B, Pikounis B, Liaw A (2003) Improved statistical methods for 
hit selection in high-throughput screening. J Biomol Screen 8: 634—647. Available: 
http://www.nebi.nlm.nih.gov/pubmed/14711389. Accessed 2011 Jul 7. 

23. Gabor EM, Alkema WBL, Jansscn DB (200 1) Quantifying the accessibility of the 
metagcnome bv random expression cloning irchniques. Environ Microbiol 6: 
879-886. Available: http://www.ncbi.nlm.nih.gov/pubmed/15305913. Ac- 
cessed 2011 Jun 22. 

24. Korecka A, de Wouters T, Gultronc A, Lapaquc N, Pettersson S, et al. (2013) 
ANGPTL4 expression induced by butyrate and rosigiitazone in human intestinal 
epithelial cells utilizes independent pathways. Am J Physiol Gastrointest liver 
Physiol 304: G1025~37. Available: htq3://www.ncbi.nlm.nih.gov/pubmed/ 
23518684. Accessed 2013 Aug 12. 

25. Gillies RJ, Didier N, Denton M (1986) Determination of cell number in 
monolayer cultures. Anal Biochem 159: 109-113. Available: http://www.ncbi. 
nlm.nih.gov/pubmed/3812988. Accessed 2011 Aug 31. 

26. Tukey JW (1977) Exploratory data analysis. Reading (MA): Addison-Wesley. 



PLOS ONE I www.plosone.org 



13 



August 2014 I Volume 9 | Issue 8 | e105598 



